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Abstract 

We compare 4d lattice results for the finite temperature phase transition in 
the SU(2)+Higgs model with 3d lattice results for the phase transition in the 
corresponding dimensionally reduced effective theory. While the large error- 
bars and the lack of a relation of the 4d lattice gauge coupling to continuum 



physics prevent rigorous conclusions, the results are nevertheless compatible. 
\Q . This provides a direct non-perturbative check of dimensional reduction in the 



present context. 



1. Due to its effect on the baryon number of the Universe, the electroweak phase tran- 
sition in the Standard Model and beyond is of considerable interest (for a review, see |l|). 
The first question to be asked is whether the phase transition is strongly enough of first 
order for cosmological consequences. It turns out that even this simple question is difficult 
to answer reliably, due to the infrared (IR) problem in finite temperature field theory. This 
has led to several lattice investigations of the problem [2-12]. 

There are two different lattice approaches available. First, there are the finite tem- 
perature 4d simulations of the bosonic SU(2)+Higgs theory [2-6]. Second, there are 3d 
simulations [7-12] of an effective theory obtained with dimensional reduction P3"|,|T4| from 



the original theory. Recently, the extrapolation to the continuum limit has been investigated 
both in 4d M and in 3d M. Hence, a comparison of the results should become possible, 
allowing to check non-perturbatively the validity of the perturbative dimensional reduction 
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used in the derivation of the effective 3d theory. We will make the comparison in a part of 
the parameter space where corrections to the 3d theory are largest and where the transition 
is strong enough to be accessible also to the 4d simulations. However, the corrections are 
still very small and hence the 4d and 3d results turn out to agree within errorbars. 

2. Let us start by briefly reviewing the main features of the two approaches. In principle, 
the more straightforward way of solving the problem is the 4d approach, in which the original 
theory is put as such on the lattice. The parameters of the lattice action are fixed from 
masses and from a suitably defined coupling constant measured at zero temperature. Then 
the lattice shape is changed so as to go finite temperature. 

Unfortunately, there are a few problems with the 4d approach. First, one cannot put 
the chiral fermions appearing e.g. in the Standard Model on lattice, and hence one can 
only study the bosonic sector of the theory. This is unacceptable, since especially the 
heavy top quark is numerically important for the electroweak phase transition. Second, it 
is not known how the renormalized lattice gauge coupling g 2 R of the 4d simulations, defined 
in terms of a static potential measured by Wilson loops, is related to continuum physics. 
Hence one does not know e.g. the MS gauge coupling g 2 (/i) at /i — m w to which some given 
4d simulation corresponds. Third, there are many different length scales in the problem 
- (2ttT)~ 1 , rriw(T) and m^ 1 (T) - - making the simulations technically demanding. In 
particular, for a large pole Higgs mass m H the inverse correlation length mJ^iT) at the 
phase transition point becomes much larger than (27rT) _1 , necessitating very large lattices. 

The 3d approach overcomes most of the direct problems of the 4d approach. Chiral 
fermions appear only in the perturbative dimensional reduction step and hence the final 
effective theory is purely bosonic. The 3d gauge coupling of the effective theory is in direct 
perturbative relation to the original zero temperature 4d gauge coupling in MS or any other 
desired scheme (in the Standard Model, it is related directly to the muon lifetime). Finally, 
dimensional reduction removes the smallest length scale (2ttT)~ 1 from the problem, making 
the simulations technically less demanding. Hence one can easily go to larger Higgs masses. 
In practice, 4d simulations have concentrated on m# = 34 GeV, whereas in 3d the continuum 
limit has been taken up to m H ps 70 GeV. 

On the other hand, the 3d approach of course relies on perturbation theory in the deriva- 
tion of the effective theory. While the perturbation theory used is free of IR-problems and 
is hence expected to be as accurate as perturbation theory at zero temperature, it might 
nevertheless be useful to check its accuracy non-perturbatively. In particular, to obtain 
the simplest possible effective theory, one neglects higher-dimensional radiatively generated 
operators with small coefficients, for instance of the type 
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>fi 



6 = c^. (1) 

In principle, there might also be non-local operators generated at high loop orders. The 
perturbative estimates made in |L4| indicate that the effects of higher-order operators should 



be on the 1% level (for other estimates, see [IB]]). 

The problems and benefits of the two approaches suggest that one should take a bosonic 
theory resembling the true original theory, and study the accuracy of dimensional reduction 
at a Higgs mass accessible to both ways of making simulations. Then one can use the 3d 
theory to include fermions and to go to higher Higgs masses. 

To investigate the accuracy of dimensional reduction, we will in this paper compare the 
results of 4d and 3d simulations for the bosonic SU(2)+Higgs model at the zero temperature 
pole Higgs mass of about 34 GeV. Hence the U(l) subgroup and the fermions of the Standard 
Model are neglected. Since the Higgs mass is small, the transition is strong and the vacuum 
expectation value v/T c = ((/>)/ T c of the Higgs field is large in the broken phase. Consequently, 
the comparison should be sensitive to higher dimensional operators of the type in eq. ([I]). 

3. The 4d results relevant for the comparison are taken from For the zero tempera- 
ture parameters, an extrapolation to the continuum limit has there been taken. The relation 
of the Higgs mass to the W mass is mn/my/ = 0.422(11), which for mw = 80.22 GeV gives 
van — 33.9(9) GeV (the number in parentheses is the error of the last shown digit). We 
shall assume m# = 34(1) GeV. The renormalized gauge coupling defined in M has the value 
9R,( m w) = 0.585(10). In principle, to get the corresponding continuum coupling one should 
calculate the static potential from which g\ is extracted also in the MS scheme. This kind 
of a calculation is not available; nevertheless, one might expect that g\{jn^) = 0.585(10) 
roughly corresponds to the MS running coupling g 2 {mw) ~ 0.585, since loop corrections 
should be of the relative magnitude g 2 /(16ir 2 ). We shall in the following consider the pos- 
sibilities that g 2 {m w ) = 0.570,0.585,0.600. 

For the properties of the phase transition, the extrapolation to the continuum limit is not 
quite perfect yet. The critical temperature has been extrapolated to be T c /m H = 2.147(40), 
but the fit is much affected by the point where the temporal extension of the lattice is 
L t = 2, relatively far from the continuum limit L t = oo. The latent heat has been measured 
to be L/T c 4 = 0.240(34) at L t = 2 and L/T c 4 = 0.28(12) at L t = 4. For the surface tension, 
values are given only for L t = 2; there two different methods give cr/T 6 3 = 0.053(5) and 
<t/T c 3 = 0.065(10) |J. For the order parameter v/T c the number 1.65(1) can be read from 
fig. 9 of , but there are large systematic errors since the relation of the lattice observable to 
continuum values is unclear (it should be noted that v is defined as v 2 oc (<p^(p), leading to a 
gauge independent but scheme dependent quantity). The 4d results used for the comparison 



3 



9 2 R ( m w) 


m H /GeV 


T c /GeV 


LlTt 




v/T c 


0.585(10) 


34(1) 


73.0(14) 

(L t = oo) 


0.28(12) 
(Lt = 4) 


0.053(5) 
(L t = 2) 


1.65(1) 

(Lt = 4) 



Table 1: The 4d lattice results for the phase transition ||||. The fit to the continuum 
limit L t = oo has been studied only for T c . To avoid confusion in the comparison with 3d 
results, the order parameter in the broken phase is denoted here by v (the definitions are 
different in 4d and 3d). 



g 2 (m w ) 


m H /GeV 


T c /GeV 


Lin 


alTl 


v/T c 


0.420 


29.1(5) 


76.8(5) 


0.200(7) 


[0.071(3)] 


1.74(3) 


0.420 


54.4(5) 


132.6(5) 


0.0294(7) 


0.0017(4) 


0.626(8) 


0.420 


64.3(5) 


151.2(5) 


0.0194(12) 


? 


0.529(20) 



Table 2: The 3d lattice results for the phase transition |J. These results represent the 
continuum limit apart from the surface tension for m# = 29.1 GeV, in which case the 
continuum limit was not investigated. 



are summarized in Table 1. 

4. The 3d results relevant for the comparison can be read from Table 10 in ||. In partic- 
ular, the second block gives the lattice values for the relevant dimensionless 3d observables, 
and the third block gives the corresponding 4d physical quantities. The results of the third 
block are reproduced in Table 2. 

Two things should be noted from Table 2, relevant for the comparison with 4d lattice 
results. First, the pole Higgs masses studied do not include 34 GeV. For the sake of the 
comparison, one should hence make an interpolation to this Higgs mass. Second, the results 
shown correspond to the MS gauge coupling g 2 (mw) = 0.420 (this value is close to that 
in the Standard Model). Since 4d simulations were made with a different gauge coupling, 
one also has to study how the results depend on g 2 . Fortunately, the scaling with g 2 can 
easily be studied since the same 3d results describe simultaneously a large class of different 
4d theories. 

We shall first explain the precise way of doing the interpolation of tuh and the scaling 
of g 2 , which was used for deriving the numerical values below. Then we also explain an 
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approximate procedure for getting a rough estimate of the numbers without having to resort 
to the precise form of the analytical formulae |14[ of dimensional reduction. 



5. The precise way of doing the interpolation of ran and change in g 2 is based on the 
dimensionless quantities x, y c , A£ 3 , £\, a 3 of the 3d theory and on their relation to 4d physics 
as explained in Sec. 11 of ||. It should be noted that the values measured for these quantities 
in 3d are independent of the value of g 2 and of the formulae for dimensional reduction. 

First, we need an interpolation of y c , A£ 3 , £ 3 , a 3 to all x. Since there are only three lattice 
points and since they are rather far from each other, we use perturbation theory to get a 
better fit. For y c which is close to zero at the critical point, we subtract from the lattice 
value the 2-loop perturbative value and then fit a parabola to the difference. For the other 
quantities, we calculate the ratio of the lattice and perturbative values and fit a parabola to 
this. With this procedure, one should get a reasonable interpolation for the critical curve 
y = y c (x) and for the values of A£ 3 , £\ along it. 

For the surface tension, there is only data point (m# = 54.4 GeV; the tree-level Higgs 
mass is then m* H = 60 GeV) in the continuum limit in Table 10 of ||, and this is pretty far 
from the point m# = 34 GeV. Hence, to get some kind of a comparison also for the surface 
tension, we include the value of a 3 at itlh = 29.1 GeV and an estimate of the surface tension 
at m* H = 70 GeV from ]10) into the fit, although these values are not extrapolations to the 
continuum limit. The point m* H = 70 GeV does not affect the result practically at all. 

Next, one needs the relations to 4d. We fix m# = 34 GeV and g 2 {m w ) = 0.570, 0.585, 
0.600. For each g 2 , we vary T and draw the corresponding curve x(rriH,T),y(mH,T) in the 
(x, y)-plane. The value of T for which this curve crosses the critical curve y = y c (x) gives the 
critical temperature. Using the 3d gauge coupling g 3 (rriH, T c ), the derivatives dx/dT, dy/dT 
at T c , the slope of the critical curve dy c /dx, and the values of A£ 3 ,£\ and a 3 at x c , one can 
then calculate v/T c ,a/Tj^ and L/T c 4 from eqs. (11.5), (11.6) and (11.7) of f§. 

With the method explained, we get the results in Table 3. We have not added any 
errors from the interpolation to the values of T C ,L/T C 4 or v/T c . The error of T c is mainly 
determined by the uncertainty in the Higgs mass, and might be slightly underestimated. For 
<t/T c 3 we give just the percentual error of the m# = 54.4 GeV continuum limit; if there were 
simulations at the pole mass m# = 34 GeV (in which case the tree-level mass parameter is 
m* H ~ 41 GeV) the errors would be considerably smaller. 

6. For a rough estimate of the numbers of Table 3, one can use a simplified procedure. 
First, just interpolate the values in Table 2 to m# = 34 GeV with the method explained 
above, giving T c = 88.4(5) GeV, L/T c 4 = 0.118(5), (x/T c 3 = 0.031(8), v/T c = 1.30(3). Then 
look from Sec. 11 of || how the values scale with g 2 = g 3 /T c w g 2 (7T c ). The critical 
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g 2 (m w ) 


m H /GeV 


T c /GeV 


l/t£ 




v/T e 


0.570 


34.0(5) 


76.7(5) 


0.19(1) 


0.048(12) 


1.43(4) 


0.585 


34.0(5) 


75.7(5) 


0.19(1) 


0.050(12) 


1.44(4) 


0.600 


34.0(5) 


74.9(5) 


0.20(1) 


0.052(13) 


1.45(4) 



Table 3: The 3d lattice results interpolated to m H = 34 GeV at different gauge couplings. 
These results represent the continuum limit apart from the surface tension (see Table 2). 



temperature is determined approximately from 

ml- -^ + cg 2 3 T 2 = 0, (2) 

where m\ is the scalar mass parameter in the dimensionally reduced 3d theory, c is a constant, 
and m 2 H is the running MS Higgs mass at a scale of order 7T C . It follows from eq. (0) that 
T c oc g% . From eqs. (11.5) and (11.6) one can see that v/T c roughly scales as g^ and 
cx/T c 3 scales as gf. To get the scaling of L/Tf, one can neglect the latter term in eq. (11.7) 
since dx/dT is non-zero only through logarithmic one-loop corrections to the dimensional 
reduction of the coupling constants. Using y = m\j g\ one then gets 

dy_ d f -m 2 H /2 + cg 2 T 2 \ _ m\ _ 2c 

where the logarithmic running of g\ and m\ with temperature was neglected and eq. (0) 
was used. Using eq. (11.7), one finally gets L/T^ oc g%- 

The values of g\ are g\ = 0.384, 0.521 for g 2 (m w ) = 0.420, 0.585, respectively. This gives 
T c = 75.9(5) GeV, L/T c 4 = 0.22(1), a/T c 3 = 0.057(15), v/T c = 1.51(4) for g 2 (m w ) = 0.585. 
These values are in rough agreement with the second row of Table 3. 

7. Let us finally compare the 4d results in Table 1 and the 3d results in Table 3. For 
the critical temperature, the 3d values are a bit higher, but compatible within errorbars for 
g 2 (mw) = 0.600. This might be an indication of the correct value of g 2 (mw)- A higher 
value of g 2 also brings the other variables closer to each other. One should also remember 
that there may be small errors e.g. in the interpolation. In addition, it should be reiterated 
that the extrapolation of the 4d critical temperature to the continuum limit might allow for 
improvement. 

The values of latent heat are clearly compatible, but the errorbars in 4d are rather large. 
The values of surface tension agree, as well, but the errorbars are again large, and neither 
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the 4d nor the 3d values represent the continuum limit at this m#. In v/T c there is some 
discrepancy which is not surprising since the definitions of v/T are different in 4d and 3d. 

We conclude that at the present level of accuracy, the results agree. This gives non- 
perturbative evidence for the applicability of dimensional reduction in the present context. 

Clearly, it would be most desirable to improve the quality of the comparison. Especially 
the extrapolation to the continuum limit should be made more precise. In the 4d case T c 
relies perhaps too much on L t = 2, and for L/T^ and cx/T c 3 a fit to L t = oo is lacking 
altogether. In the 3d case the fit to the continuum limit for cr/T c 3 at m# = 54.4 GeV 
(m* H = 60 GeV) contains large errors since the transition is already rather weak there; 
for rriH = 29.1 GeV {m* H = 35 GeV) the continuum limit is much easier but was not 
studied in ||. Of course, it would also be most important to get some direct evidence 
on the value of g 2 (mw) to which the 4d simulations correspond. A comparison of finite 
temperature correlation lengths, which seems not to be available at the moment, would be 
welcome. Finally, all the comparisons would be more straightforward and the errors related 
to interpolation removed if lattice simulations were made with the same pole Higgs masses. 
Once the errors have been reduced, the effects of higher-order operators can be reliably 
estimated. 
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